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, We study the unpolarised and polarised hadro-production of charmonium in non-relativistic QCD (NRQCD) 

■ at low transverse momentum, including sufficiently higher orders in the relative velocity, v, so as to study the 
^\ ' ratio of xa and Xc2 production rates. 
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1. Introduction 

Recent progress in the understanding of cross 
sections for production of heavy quarkonium res- 
onances has come through the NRQCD reformu- 
lation of this problem A Factorisation ap- 
proach based on NRQCD developed by Bodwin, 
Braaten and Lepage |^ enables the factorisation 
of the production cross sections for a quarkonium 
H (with momentum P) into a perturbative part 
and a non perturbative part — 



da 



where $ is a flux factor and 11 (i?) denotes the 
hadronic projection operator. The fcrmion bi- 
linear operators ICi are built out of heavy quark 
fields sandwiching colour and spin matrices and 
the covariant derivative D. The labels i, j include 
the colour index, spin S, orbital angular momen- 
tum L (coupling the N covariant derivatives), the 
total angular momentum J and helicity Jz . 

The coefficient function Cij is computable in 
perturbative QCD and has an expansion in the 
strong coupling as{m) (where m is the mass of 
heavy quark), whereas the matrix elements are 
non perturbative. However, in NRQCD, these 
matrix elements scale as powers of v. Hence the 
resulting cross section is an expansion in powers 
of as(rn) and v. Often, higher orders in v in- 
volves previously neglected colour-octet states of 
the heavy quark pairs. For charmonium states, 
a numerical coincidence, ~ q;s(to^), makes 
the higher order terms in v important and the 
double expansion more complicated. The above 
non perturbative matrix elements can be reduced 
to the diagonal form, 0^(^'^+^Lj') (where a de- 
notes the colour singlet or octet state) and oS- 



diagonal form, 75^(2^^+1 Ly,2^?+iLy') g. These 
matrix elements scale as v"^ where d = 3+N+N'+ 
2{Ed + 2Md), Ed and Md are the number of colour 
electric and magnetic transitions. This formal- 
ism has been successfully applied to large trans- 
verse momentum processes . Inclusive produc- 
tion cross sections for charmonium at low ener- 
gies, dominated by low transverse momenta, also 
seem to have a good phenomenological descrip- 
tion in terms of this approach |^-^. The spin 
asymmetries have also been computed both for 
the low transverse momentum processes H] and 
for those with high transverse momenta 

It was argued in |^ that a better understand- 
ing of such cross sections, and asymmetries, can 
be obtained if the higher order terms in v and as 
are used. This follows from the fact that the total 
inclusive J/ip cross sections arise either from di- 
rect J/ip production (which starts at order asv"^) 
or through radiative decays of XJ states, xo a-nd 
X2 are first produced at order asV^, whereas xi; 
which has the largest branching fraction into J/V', 
is produced only at order UsV^. Further phe- 
nomenological problem is to explain the Xi/X2 ra- 
tio observed in hadro production [^,0 . A better 
understanding of these cross sections requires the 
NRQCD expansion up to order asV^. The unpo- 
larised and polarisation cross sections are defined 
as 

a = ^a{h,h'), Act = ^ ;i/i'ct(/i, /i'), (2) 



hh' 



hh' 



where h and h' are the helicities of the beam and 
target respectively, and a{h, h') denotes the cross 
section for fixed initial helicities. The difference 
between the polarised and unpolarised case is in 
the coefficient functions denoted by Cij for the 
polarised case; the set of non perturbative matrix 
element are the same. We construct the coeffi- 
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cient functions using the "threshold expansion" 
method ||Tl| ] and enumerate the non perturba- 
tive matrix elements using the spherical tensor 
method described in 11121. 



2. Cross Sections 

To lowest order in as the contributing parton 
level cross section are qq QQ and gg — > QQ. 
The hadron level cross section is obtained by mul- 
tiplying appropriate parton luminosities — 



a{xi) b{x2), 



(3) 



where a, b runs over quark qj, antiquark g/ or 
gluon g densities depending on the subprocess 
cross sections. For polarised cross sections, Act, 
the corresponding polarised luminosities ACab is 
obtained by replacing a, by polarised parton den- 
sities Aa. Data indicates that £qq <C Cgg for 
> 20 GeV and \ACqq\ < Cgg. Consequently, 
the qq channel may be neglected for double po- 
larised asymmetries to good precision. 

The squared matrix element for the gg process 
is technically more complicated. The difference 
between the unpolarised and polarised cases 
|l3|| lies solely in the flipped sign of the J = 2 
part, which arise in the polarisation sum of initial 
state gluons. The subprocess cross section for the 
production of a charmonium H can be written as 



4m? 



S{s 



4m^ 
5 



18 



48 



16 



(4) 



where d runs over the various matrix elements 
that contribute to the charmonium H at order 
v'^. The subscript S, D and F denotes the colour 
singlet, colour octet symmetric and antisymmet- 
ric parts respectively. For the polarised case 
is replaced by Aa^ and the combination of non 
perturbative matrix elements by . For the 
various charmonium states J/iJj and XJ the com- 
bination of non perturbative matrix elements are 
listed bellow. The changes for the polarised case 
is mentioned appropriately. 



Direct J/tp Production 



2m? 
1 

'2^ 
1 



35 



15mS 

I .-'/'/' ,'3 pi 3i 



(5) 



-jivr{^pl^pl) + 2V'r{'Ph'Pi) 



1 



2m^ 
Xo Production 

3 



ef{5) 



ef{7) 



2m4 
7^/5 

4\/3m' 



On'Po'), 



ef{9) 



1 



8m 
U9V7 



245 



Of (3p3) 



Pf"(3p„l,3p5) 



10V3 



5m' 



-ore Pi), 



(6) 



2m^ 
1 

2^ 
1 



30f(3pj) + -0f eP^i) 



18m6 

Xi Production 
1 



ef{9) = 

01^(9) = 



ef{9) 



2m4 
1 

2m^ 
1 

'2^ 
1 



30f (3Poi) + -Of (3p,i) 



30f (3pJ) + -Of ('^^2^) 



(7) 



5,^A/1/^^o^ 1 



Gm'^ 



3m6 



In eqs.(^,^,^ the coefficient of J = 2 matrix ele- 
ments changes sign for polarised cases. 

Xi is produced first at order v^. The large 
branching ratio for the decay xi ~^ J /''P makes 
this a phenomenologically important term, and is 
the main motivation for this work. 
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X2 Production 

er(7) . 



15m' 



Of (3p,l), 



Of (9) 



eF(9) 



-Of (3p„i) + 



262 



ore Pi) 



2V7 
1 



2m^ ^ 
1 

2^ 
1 



Of (1^0) 



1 



^of r5f) + ^of (^i^f) 
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In the combination 0f (9), the coefficient of 
Of (■^Pf) is different for the polarised case, both 
in magnitude as well as sign. This is because 
the matrix element arises also from sources other 
than the polarisation sum that contribute to the 
J = 2 part. To obtain the polarised expression 
replace the coefficient 262 by —238. The other 
J = 2 terms change sign as expected. The J = 3 
term in the F colour amplitude is also derived 
from the J — 2 part of the polarisation sum, and 
hence ffips sign compared to the unpolarised case. 

3. Discussion 

In spite of the large number of unknown non 
perturbative matrix elements in the final re- 
sults, it is possible to make several quantitative 
and qualitative comments about the polarisation 
asymmetries by making use of heavy quark spin 
symmetry and scaling arguments developed in 
13 1 . This scaling argument allows us to make 
rough estimates. Neglecting possible logarithms 
of m and dimensional argument can be used to 
write 



{lC^Ii{H)lC] 



(9) 



where Dij is the mass dimension of the operator, 
d is the velocity scaling exponent in NRQCD, A 
is the cutoff scale below which NRQCD is defined 
(A ~ m), and Yij and are dimensionless num- 
bers. Rh contains the irreducible minimum non 
perturbative information. 
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Figure 1. The ratio (j{xi) I o{x2) measured inpN 
and ttN has been plotted as a function of \/S. 
The data has be compiled in Preliminary re- 
sult presented by E771 at the Moriond conference 
is 0.45 ± 0.2 |1^ and in agreement with our esti- 
mates. 



Assuming the constancy of Rh fl^ and us- 
ing heavy-quark symmetry, we find that the non 
perturbative matrix elements contribute approx- 
imately 12R^m^v^ to the xi cross section and 
about (1 + + 12v*)R^m'^v^ to the X2 cross sec- 



tions. Then we expect 



1 



12^4 



0.45, 



(10) 



independent of \fS. This estimate is in rea- 
sonable agreement with the measured values in 
proton-nucleon collisions (see fig. 1) — 0.34 ±0.16 
atVS^ 38.8 GeV [0 and 0.24 ± 0.28 at VS = 
19.4 ||l^. The measurements are also compat- 
ible with a lack of ^/S dependence. Estimate 
of 0{al) effects in Q was used to show that 
the xi/X2 ratio could be about 0.3. In NRQCD 
this ratio cannot depend on the beam hadron. 
It turns out that the estimate in eq. (jl^) is not 
very far from the recently measured value in pion- 
nucleon collisions — 0.57 ± 0.19 at = 31.1 
GeV However, the experimental situation 

certainly needs clarification. 
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A straightforward application of the NRQCD 
scaUng laws would lead us to the conclusion that 
the asymmetries for pp — > xo.2 are given by 



AXa 
pp 



— A^^ 
pp 



AC, 
C 



(11) 



aa 



where we have neglected the contribution of the 
qq channel. The asymmetry for xi production is 



AXl 

pp 



(9) 



0^)^(9) 



6^,^(9) 



Co 



(12) 



The ratio of matrix elements can be estimated us- 
ing heavy quark spin symmetry and the scaling 
relations in eq. (^). 0^'(9) vanishes in this ap- 
proximation and the terms Q^''p{9) come with op- 
posite signs. The numerator is positive but small 
and we expect — 

AC 

Q.2=P^. (13) 



4X1 
pp 



The qq channel remains negligible even at ^/S — 
500 GeV. 

The J 1^1) asymmetry seems to be enormously 
complicated because of the radiative decays of the 
X states. However, a major simplification occurs 
because of the near vanishing asymmetry in direct 
J 1^1) production. Thus the asymmetry comes en- 
tirely from the 20-40% of the cross section due to 
X decays. Taking into account the ratios of the 
production cross sections of x and the branch- 
ing fractions for their decays into J/V', we find 
that the xi and X2 states contribute equally to 
J/'tp. Hence the J/ if) polarisation asymmetry is 
expected to be approximately 

AC, 



-(0.15 ±0.05)- 



-93 



c. 



(14) 



We summarise the predictions made on the basis 
of the NRQCD scaling in eq. (^) and the assump- 
tion of Ku depends only on the hadron H — 



■t^pp ~ J^pp ^pp 



_^X2 
PP 



A£„ 



(15) 



In conclusion, low energy double polarised 
asymmetries are a good test-bed for understand- 
ing the origin of all observed systematics in fixed 
target hadro-production of charmonium. The 
high order computations presented here provides 
a set of processes which can be used to test as- 
pects of NRQCD factorisation and scaling. 
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